Heat conduction in the disordered Fermi-Pasta-Ulam chain 



00 

o 
o 

(N 

in 



I 

CZ3 



I 

O 

o 



> 

\o 
o 

o 

oo 

o 



Abhishek Dhar^'B and Keiji Saito^^^ 

'Raman Research Institute, Bangalore 560080 
'^Graduate School of Science, University of Tokyo, 113-0033, Japan 
^CREST, Japan Science and Technology (JST), Saitama, 332-0012, Japan 
(Dated: June 25, 2008) 

We address the question of the effect of disorder on heat conduction in an anharmonic chain with 
interactions given by the Fermi-Pasta-Ulam (FPU) potentiaL In contrast to the conclusions of an 
earlier paper [Phys. Rev. Lett. 86, 63 (2001)] which found that disorder could induce a finite 
thermal conductivity at low temperatures, we find no evidence of a finite temperature transition 
in conducting properties. Instead, we find that at low temperatures, small system size transport 
properties are dominated by disorder but the asymptotic system size dependence of current is given 
by the usual FPU result J ~ . We also present new interesting results on the binary-mass 

ordered FPU chain. 

PACS numbers: 



It is now generally believed that heat conduction in 
one-dimensional (II?) systems is anomalous jll, In the 
absence of an external pinning potential, as is the case in 
most realistic situations, one finds that Fourier's law is 
not satisfied. One of the predictions from Fourier's law 
is the scaling form of heat current J with system size N 
for a system with a fixed applied temperature difference. 
From Fourier's law one gets J ^ . The conclusion, 
from a large number of studies of ID momentum con- 
serving systems seems to be: 
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The main results on the exponent a can be summarized 
as follows: (1) for a pure harmonic chain J ~ iV" Q, 
(2) for a disordered harmonic chain a depends on the 
spectral properties of the bath and on boundary condi- 
tions 0, H, 01 1 and (3) for a nonlinear system without 
disorder a seems to be independent of properties of heat 
baths and the results from the most recent simulations 
indicate a universal value of a = 1/3 [1, 0] ■ 

Transport in systems with both disorder and interac- 
tions has recently attracted a lot of interest both theoret- 
ically [1, 0, [l^, 11 , l2| and experimentally [l^l . The main 
interest is to understand the transition from an insulating 
state governed by the physics of Anderson localization to 
a conducting state as one increases interactions. In the 
context of oscillator chains we note that the physics of the 
disordered harmonic chain is dominated by localization 
physics which has its strongest effect in \D systems. The 
question of the effect of anharmonicity on localization 
was recently addressed for a system where the harmonic 
part of the Hamiltonian included an external pinning po- 
tential and the anharmonicity was a quartic onsite poten- 
tial [ll[. In this case, in the absence of the anharmonic 
term, J ~ e~'^^ . Surprisingly it was found that adding 
a small amount of anharmonicity leads to a conducting 
(Fourier- like) behaviour with a power law decay J ^ XjN 
and no transition to the insulating state was found on de- 



creasing the anharmonicity. An important feature seen 
was that with decreasing anharmonicity one had to go to 
larger system sizes to see the true asymptotic behaviour 
of the current. 

In the present letter we investigate the same question, 
namely that of the effect of interactions (phonon-phonon) 
on localization, but in the absence of any external pin- 
ning potential. We study the mass disordered FPU model 
with interactions put in through a quartic inter-particle 
potential. In the absence of pinning, for the harmonic 
case, low frequency modes with w l/N^/"^ remain ex- 
tended and this gives rise to a power law dependence of 
J on Af 0, [i, [i]. On the other hand, for the pure FPU 
chain also, low frequency modes are believed to play an 
important role in transport and give rise to anomalous 
transport. An earlier study by Li et al 0] on this model 
concluded that this model showed a transition, from a 
Fourier like scaling J ^ 1/A^ at low temperatures, to 
a pure FPU Uke behaviour with J - 1/A^"-^^ at high 
temperatures. Our study suggests that this conclusion 
may not be correct. We do not find any evidence of 
a finite temperature transition. Instead we find that a 
small amount of anharmonicity leads to the same system 
size dependence as seen in the pure system. We discuss 
possible sources of error in the conclusions of Li et al. 
We also present new results on the ordered binary mass 
FPU chain including nontrivial scaling properties of the 
system size dependence of current. 

Model. — We consider the following FPU Hamiltonian: 
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where {xi , pi } denotes the position and momenta of par- 
ticles and we use fixed boundary conditions xq — xn+i = 
0. The interparticle harmonic spring constant has been 
set to one and v denotes the strength of the quartic in- 
teraction. We consider a binary random alloy and set 
the masses of half of the particles, at randomly chosen 



2 



sites, to mi and the rest to m2- The particles at the two 
ends of the chain are connected to stochastic white noise 
heat baths at different temperatures. The equations of 
motion of the chain are then given by: 

mixi = -{2xi - xi-i - xi+i) 

- v[{xi - xi^i)^ + {xi - xi+i)^] - -fixi + r]i , (3) 

with T]i = ?7l<5;,i + r]R^i,N, li = + Si^n), and where 

the noise terms satisfy the fluctuation dissipation rela- 
tions {vL{t)vL{t')) = 2-fkBTL5{t - t'), (??fl(t)?7fl(t')) = 
2'ykBTji6{t — t'), ks being Boltzmann's constant. The 
heat current is given by J = ^i{fi,i-iXi) /{N — 1) where 
fi,i-i is the force exerted by the (/ — l)th particle on the 
Ith particle and (...) denotes a steady state average. We 
will denote by [J] an average over disorder. As noted 



in [li], Eqs. ([3]) are invariant under the transformation 
Tl,r sTl^r, {xi} ^ {s^/'^xi} and v ^ vjs. This im- 
plies the scaling relation J{sTl, sTr, v) = sJ{Tl, Tr, sv) 
and thus the effect of changing v can be equivalently 
studied by changing T^, T^. In the present study we will 
fix temperatures and consider the effect of changing i/. 

Let us consider first the harmonic case i/ = 0. In this 
case it is known from detailed numerical work and an- 
alytic arguments that the exponent a depends on the 
properties of the bath and on boundary conditions. For 
white noise baths one finds a — —1/2 for fixed bound- 
aries and a = 1/2 for free boundaries. In the presence 
of anharmonicity it is expected, and indeed we have veri- 
fied in simulations, that a does not depend on boundary 
conditions. Here we use fixed boundary conditions only. 

Before presenting the results of simulations for the 
binary-mass disordered anharmonic chain it is important 
to know the value of a for the binary ordered chain. Let 
us thus discuss this first. This case was earlier studied 
in 0] where it was found that the temperature profile 
showed oscillations with an amplitude that decreased as 
N~^/'^. Let us denote the mass ratio mi/m2 = A. For 
A — 1^ the simulations in Q gave strong evidence for 
an exponent a = 1/3. However for the value A — 2.62 a 
clear convergence could not be obtained. Here we will ar- 
gue that the exponent remains unchanged from the A = 1 
value. In Fig. ([1]) we present simulation results for the iV- 
dependence of the current J in the binary-mass ordered 
chain for different values of the parameter A, all corre- 
sponding to the same average mass (toi -I- m2)/2 = 1. 
Remarkably we find that at large enough system sizes 
the actual values of the currents for different A tend to 
converge to the same value as the A = 1 value. Thus 
clearly the exponent a remains unchanged for any value 
of A. However for a large mass ratio one has to go to 
large system sizes to see the true exponent. A similar ef- 
fect was seen in Refn.Ql for the binary hard-particle gas. 
In our simulations we used the velocity- Verlet algorithm 
with time steps dt = 0.005 [l^. For small system sizes 
we used O(IO^) — O(IO^) steps for relaxation and same 
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FIG. 1: Plot of the heat current J versus system size in the 
binary mass ordered chain for different values of the mass 
ratio A = 1.0, 1.1, 1.22 and 1.5. 



number of steps for averaging, while for larger systems, 
up to O(IO^) steps were used. In all our simulations we 
used Tl = 1.25, Tij = 0.75 and 7 = 1.0. 

For the disordered anharmonic case, we wish to study 
the cases with weak and strong anharmonicity and see if 
there is a transition in the value of a. For small v we 
expect the system's behaviour to be close to a harmonic 
one, and so one would have to go to very large system 
sizes to see the effect of anharmonicity and the correct ex- 
ponent. It then becomes necessary to try and understand 
the data using some sort of a scaling analysis. Let us first 
do this for the ordered case. We fix the value of ^4 = 1.5 
and look at the A^-dependence of the current for different 
values of v. The results are shown in Fig. ((2^). For small 
system sizes, we find a flat region which is expected since 
for system sizes smaller than the phonon-phonon scat- 
tering length scale we expect the system to behave as a 
harmonic chain. The scattering length should be larger 
for smaller f and this can be seen in the plot. At large 
enough system sizes all curves tend to show the same de- 
cay coefficient a = 2/3. To see this clearly we scale the 
system size by a length factor Fig. ^p) shows a 

nice collapse of the data and the value of the exponent is 
confirmed. We find empirically that the i^-dependence of 
the length parameter is given by £{i') = l/tanh(2i/). We 
note the interesting and somewhat surprising point that 
for any given system size, the value of the current satu- 
rates as we keep increasing i/. In Fig. ([3^) we show the 
temperature profiles for different system size for ^' = 0.1. 
As noted earlier in [3] we see the large oscillations in 
temperature. An interesting general feature of temper- 
ature profiles in FPU chains is the following. A coarse- 
grained temperature profile obtained by averaging over 
many particles would be smooth and monotonic. How- 
ever the temperature gradient is non-monotonic and this 
appears to be true even for small temperature differences. 
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FIG. 2: Plot of the heat current J versus system size for the 
binary mass ordered chain, for different values of v and with 
the mass ratio A = 1.5 (a). Fig (b) shows the same data 
plotted with a scaled a;-axis. 



This implies that it may not even be possible to write a 
phenomenological relation such as J = —kn(T)\7T. 

Finally we now give the results for the disordered an- 
harmonic case. We take averages over 50 — 100 samples 
N < 1024, 10 samples for N = 1024 - 16384, and 2 sam- 
ples for N = 32768 and 65536. In Fig. (g]) we plot the 
results of simulations for [J] for v — 0.004 and 0.02. Also 
we show the result for = 0.0. For small values of we 
see that, at small system sizes the current value is close 
to the i/ = value. As expected one has to go to large 
system sizes to see the effect of the weak anharmonicity. 
At sufficiently large N we find the same system size de- 
pendence of [J] as obtained for the ordered FPU chain, 
namely that given by a = 1/3. In fact by scaling the 
current by appropriate factors we find that the data for 
the disordered case can be made to collapse on to the 
binary-mass ordered case. This is shown in Fig. ([5]) (for 
v = 0.02, 0.1, 2.0). Thus our results show that the asymp- 
totic power law dependence of the current is always domi- 
nated by anharmonicity while disorder only decreases the 
overall conductance of a sample. In Fig. ([SJs) we plot the 
temperature profile for v = 0.1 and find that the asymp- 
totic profile is similar to the ordered case. We now throw 
some light on the reasons which led to the erroneous con- 
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FIG. 3: Plot of the temperature profiles in the (a) ordered 
and (b) disordered lattices for v = 0.1 and for different system 
sizes. 



elusions in [8|, of a transition in conducting properties 
at low temperatures (or equivalently small anharmonic- 
ity). Consider the data for [J]N plotted in Fig. (|4]) for 
v = 0.004. We see that at around N ~ 1000 - 2000 the 
data seems to flatten and if one had just looked at data in 
this range, as was done by [1], one would conclude that 
Fourier's law is valid. However the behaviour changes 
drastically when one looks at larger system sizes and one 
again gets the usual FPU behaviour. To verify that this 
is indeed what happens for the particular case studied by 
Li et al we have repeated simulations with their set of pa- 
rameters but for much larger system sizes and the results 
are shown in the inset of Fig. This case corresponds 
to a much smaller value of v and so it is expected that 
it will follow the i' — curve till very long length scales 
and this is clearly seen. However at around N = 16384 
we see a tendency for the curve to turn up and we expect 
that the same asymptotic behaviour to eventually show 
up. While a transition cannot be ruled out at even lower 
temperatures, this seems unlikely. Also, if there is such 
a transition, it should probably be to a disordered phase 

with [J] - l/A^3/2_ 

Temperature dependence of conductivity. — The scaling 
property of the current, mentioned earlier [after Eq. ([3])], 
implies that the thermal conductivity has the form k ~ 
k{i'T). For small anharmonicity {v << 1), our earlier 
results imply that at large system sizes k ~ (Ni^y^^ 
and from the scaling property this immediately gives 



4 



50 
20 

^ 5 

2 
1 



0.01 

IJ]N 

0.001 



X WN 
+ iy = d.O 
□ , NH , 



10^ 10' 10"" 



* i 



" X = 0.02 
V = 0.004 
□ u = 0.0 , 

' I ' ' I I I I I 



—I I I i_ 



10^ 



lO'' 

N 



10* 



FIG. 4: Plot of heat current versus system size, for the dis- 
ordered anharmonic chain, for different values of u. The 
data in the inset, corresponds to parameters {Tl,Tr) — 
(0.001, 0.0005) with Gaussian white noise bath for u — 1 
(WN) and v = 0, and Nose-Hoover bath (NH) for z/ = 1. 
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FIG. 5: Plot of scaled heat current versus scaled system size, 
for the binary-mass ordered (J) and disordered ([J]s) anhar- 
monic chains, for different values of u. 



K ^ T^/^ at low temperatures. However at small system 
sizes [A^ << ^{v)], we expect the system to behave like 
a harmonic system with k ^ T^. At high temperatures 
the conductivity will saturate to a constant value. 

Discussion. — Our main conclusion is that there is no 
change, in the asymptotic power law dependence of the 
current on system size, on decreasing the temperature in 
the disordered FPU problem. At low temperatures one 
has to go to much larger system sizes to see the true ex- 
ponent, whose value (a = 1/3) is the same as that for 
the ordered FPU chain. We also find several interesting 
new results for the binary-mass ordered FPU chain: (i) 
the exponent a is independent of the mass ratio A and 
is the same as the A — 1 value, (ii) the data for different 



values of v can be collapsed by scaling the system size by 
a ;y-dependent length scale. Also we make the interest- 
ing and somewhat surprising observation that for a finite 
FPU chain, J — > const 7^ as z^ — ^ oo. Experiments 
measuring heat conduction in quasi- ID systems such as 
nanowires, nanotubes [l^ are now being done. The ef- 
fect of isotopic disorder has also been measured [l^ . Our 
prediction is that while disorder will lower the current in 
a wire, the system-size dependence of [J] is unaffected. 
Experimentally, the temperature dependence of the ther- 
mal conductivity may be easier to measure and one can 
verify if this is unaffected by disorder. 
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